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Unimodular gravity is classically equivalent to General Relativity. This equivalence extends to 
actions which are functions of the curvature scalar. At the quantum level, the dynamics could differ. 
Most importantly, the cosmological constant is not a coupling in the unimodular action, providing 
a new vantage point from which to address the cosmological constant fine-tuning problem. Here, a 
quantum theory based on the asymptotic safety scenario is studied, and evidence for an interacting 
fixed point in unimodular f{R) gravity is found. We study the fixed point and its properties, 
and also discuss the compatibility of unimodular asymptotic safety with dynamical matter, finding 
evidence for its compatibility with the matter degrees of freedom of the Standard Model. 


I. INTRODUCTION 

There are only very few observations from which one could expect to learn something about the deep structure 
of spacetime, described by a model of quantum gravity. One of those is the observed accelerated expansion of the 
universe, which can be modelled by including a nonzero cosmological constant in the Einstein-Hilbert action. This 
entails the cosmological constant problem. Here, we will focus on one aspect of this problem, namely the question why 
quantum vacuum fluctuations do not seem to gravity, i.e., why the cosmological constant exhibits a severe fine-tuning 
problem. As its mass-dimensionality is 2, one would expect quantum fluctuations to drive it to be of order one in units 
dehned by the square of the physical mass scale of the theory, which is the Planck scale. In terms of Renormalization 
Group (RG) trajectories, the tiny value inferred from observations implies that a particularly fine-tuned trajectory 
has to be picked. Of course every so-called relevant coupling in a quantum theory corresponds to one free parameter 
that can only be hxed by comparison with experiment. Thus one always has to pick a particular RG trajectory in 
order for the model to reproduce observations. On the other hand, relevant couplings with only, e.g., a logarithmic 
running imply that if one picks a trajectory nearby, the measured value of the coupling will only change slightly. For 
the cosmological constant, this statement is not true, i.e., for a reason which as yet has no dynamical explanation, 
our universe just ’’happens” to live on a highly fine tuned choice of trajectory. 

As already proposed by Weinberg [T], a ’’degravitation” of the cosmological constant is possible by changing its 
status from a coupling in the action to a constant of integration that arises at the level of the equations of motion. 
The second is a ’’classical” quantity in the sense that it is not affected by quantum fluctuations. In unimodular 
gravity [ 5 ], the metric is conceived as a symmetric tensor with fixed determinant -y/—g = e [ 3 ]- This implies that no 
operator of the form y/—g exists, as the volume is just a fixed number, and the cosmological constant is removed from 
the space of couplings, the theory space. Once a Renormalization Group trajectory in this reduced space has been 
picked, the effective equations of motion can be calculated from the full effective action - the infrared endpoint of the 
trajectory - and the cosmological constant will then make its appearence as a constant of integration. In this way, 
quantum vacuum fluctuations do not affect the value of the cosmological constant. It is thus of interest to investigate 
a quantum theory of unimodular gravity. Unimodular gravity in both its quantum and classical form has sparked 
considerable interest since it was originally proposed [3HS]. 

A second motivation to consider unimodular quantum gravity lies in the fact it will most probably differ from 
the non-unimodular version of quantum gravity. This inequivalence arises, as imposing the unimodularity condition 
alters the spectrum of fluctuations of the theory. In more detail, deriving the full metric propagator by taking the 
second variation of the action yields different results when y/—g = e is imposed, than if the metric determinant is 
allowed to fluctuate. In particular, fluctuations of the conformal mode, which yield an instability of the path-integral 
in the Euclidean case when starting from the Einstein-Hilbert action, are absent in unimodular gravity. This already 
suggests that although classically equivalent [^, the quantum version of unimodular and ’’standard” gravity could 
differ. The absence of the conformal instability even suggests that the unimodular quantum theory could have better 
properties. 

As the search for the ’’quantum theory of gravity describing our universe” is still ongoing, an exploration of different 
models for quantum gravity models is clearly of interest, both from a theoretical as well as from a phenomenological 
point of view. Here we will focus on exploring models of asymptotically safe quantum gravity. As we will point out 
in Sec. |Vj it might be possible to distinguish between different versions of asymptotically safe gravity experimentally. 

Finally, in order to better understand the structure of Renormalization Group (RG) flows in gravity it is of interest 
to consider settings with fewer propagating degrees of freedom in the path-integral. Here, we should clarify that the 
physically propagating degrees of freedom in both settings, unimodular vs. ’’full” gravity, agree in a perturbative 
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expansion around a flat background, i.e., there is a massless spin-2 graviton |5]. On the other hand, the configurations 
that enter the path integral are different metric configurations, and, e.g., the configuration space of the conformal 
mode is also summed over in the case of ’’full” gravity. In order to shed light on the physical mechanism of asymptotic 
safety, it is helpful to consider settings where some of the modes in the path-integral are removed. 

In this paper, we will consider unimodular asymptotic safety, and investigate truncated Renormalization Group 
flows based on an f{R) action. We will discuss the classical equivalence of ’’full” gravity, which we will call Einstein 
gravity, with unimodular gravity based on an /(i?) action in Sec. IIA We will then focus on the quantum theory, 
and review the asymptotic safety scenario and the functional RG in Sec. m where we will also present all technical 
details of our calculation. In Sec. IV we will present the flow equation for f{R) and discuss a fixed point and its 
properties. We will make a first step toward phenomenology in studying the effect of dynamical matter in Sec. |Vj 
and finally conclude in Sec. |V^ 


II. RELATION BETWEEN UNIMODULAR GRAVITY AND EINSTEIN GRAVITY 
A. Glassical equivalence of f{R) gravity and unimodular f{R) 

Before we embark on an analysis of quantum gravity, let us clarify the classical relation between unimodular f{R) 
gravity and f{R) gravity with a full metric, see, e.g., [10] for a review. Here we will focus on the Lorentzian case, and 
then switch to a Euclidean setting for the analysis of the quantum theory. We focus on actions given by a function of 
the curvature scalar, f{R), with /(O) = 0. We introduce the Newton coupling Gy and the cosmological constant A 
in the action 


S = I d^x^g (^f{R) + ^^A + 


Cr. 


The corresponding equations of motion are given by 


- + f[R)R^u - D^D^fiR) + g^.D^fiR) + ^ 

Z lOTTLrjV 


R9p,v — „ Rfivi 


IfiTTGAr""^" - 2 


where the energy-momentum tensor is given by 


The Bianchi-identities 


^ _ 2 6L^ 


( R^,v - ^9 )i.vR ) = 0 , 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


will now play a crucial role: By taking the covariant derivative of Eq. ([^, we deduce the conservation law for the 
energy-momentum tensor by imposing the Bianchi-identities. 

To obtain the unimodular equations of motion for the action 


Su = j d'^xeU{R)+C^), 

we have to consider tracefree variations g^'^Sg^i, = 0 and obtain 

riR)R^, - D^D.f'iR) + ^g^^D^f'iR) - ^g^.RfiR) = ^ • 

Crucially, the covariant derivative of the Ihs of Eq. does not vanish when the Bianchi-identities are used. Instead, 
we can impose conservation of the energy-momentum tensor^, and thereby derive a nontrivial identity, namely 


( 5 ) 


( 6 ) 


D, 




( 7 ) 


^ As discussed, e.g., in [3, it is not clear whether this requirement can be preserved in a quantum field theory setting for the matter 
degrees of freedom. This could potentially lead to a situation where the low-energy effective equations of motion allow us to distinguish 
between General Relativity and unimodular gravity. 
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where we have used Eq. ^ and imposed Eq. Q. This identity allows us to identify l/2/(i?) with —jf'{R)R — 
jD^f'{R) + \T^, up to a constant of integration, which we choose to call A. Inserting this identity into Eq. ([^ 

we obtain Eq. (il.i.e., classically unimodular f{R) gravity cannot be distinguished from standard f{R) gravity. Note 
that this statement depends on the postulate of energy-momentum conservation in the unimodular case. It is a priori 
clear that the two theories can only be classically equivalent, if one additional condition is imposed in the unimodular 
case: Since the equations of motion of unimodular gravity are obtained by removing the trace from the standard 
equations of motion, they contain precisely one condition less. 


III. UNIMODULAR QUANTUM GRAVITY 
A. Asymptotic safety 

From now on we will focus on Euclidean quantum gravity, as this allows for a straightforward application of RG 
tools. To arrive at a unimodular quantum theory of gravity, we will invoke the asymptotic safety conjecture for gravity 
m Interestingly, unimodularity plays a role in a number of other approaches to quantum gravity, e.g., within causal 
set quantum gravity, where a discrete version of unimodularity would be implemented by performing the path-sum 
over all causets with a fixed number of elements nanni. Further, Causal Dynamical Triangulations is based on 
a setting where the number of simplices usually is held fixed for the simulations, i.e., the cosmological constant is 
removed from the space of couplings, see, e.g., El- 

An asymptotically safe quantum theory of gravity is valid at arbitrarily high momenta, i.e., beyond the regime of 
validity of effective field theory and at the same time remains predictive, i.e., only comes with a finite number 

of free parameters. It is the second requirement that breaks down in a perturbative quantization of gravity when 
extended to arbitrarily high momenta [IIHIS]- Then, an infinite number of counterterms needs to be introduced and 
there is no mechanism to predict the corresponding free parameters. In asymptotic safety, the Renormalization Group 
flow approaches an interacting fixed point at high momenta. Specifically, we refer to the dimensionless couplings here, 
which can be obtained from the dimensionful ones by an appropriate rescaling with the RG momentum scale k. If these 
couplings approach a fixed point, the theory becomes scale free and can be extended to arbitrarily high momentum 
scales. The RG flow lives in an infinite dimensional space of all couplings, the theory space, as quantum fluctuations 
generically generate all operators that are compatible with the symmetries. In this infinite-dimensional space one 
must then investigate whether predictivity can be obtained, i.e., whether the model has only a finite number of free 
parameters. This is ensured, if the interacting fixed point comes with a finite number of relevant, i.e., ultraviolet 
(UV) attractive directions. The (ir)relevance of a coupling gi{k) determines its scale-dependence in the vicinity of a 
fixed point. This can be obtained after linearising the RG flow: 

/ h \ 

. ( 8 ) 

Herein, denotes the fixed-point values of the coupling gi. are the eigenvectors of the stability matrix M.ij = 
{dj3gjdgj)\g^=g^,, uud —0i its eigenvalues. C/ are constants of integration. If Oj < 0, then Cj = 0 is required in 
order for the couplings to approach the fixed point in the UV limit, where A: —)■ oo. On the other hand, relevant 
directions are those with 9j > 0. These approach the fixed point automatically, imposing no requirement on the 
corresponding parameter Cj. They have instead to be determined experimentally, kg is an arbitrary reference scale, 
which can be taken as the scale at which the value of the couplings is measured. At a Gaufiian fixed point, such as that 
underlying asymptotic freedom in Yang-Mills theory, only the couplings of positive and vanishing mass dimensionality 
can be relevant. At an interacting fixed point, quantum fluctuations shift the critical exponents away from the mass 
dimensionality by an anomalous dimension. Additionally, the (ir)relevant directions are typically no longer given by 
the operators defined at the Gaufiian fixed point, but in fact correspond to mixtures of these. When the anomalous 
dimensions remain finite, this suggests that only a finite number of relevant directions exist, and that these can be 
found among appropriate combinations of the couplings with the largest mass dimensionalities. 

It remains for us to determine whether an interacting fixed point exists in unimodular gravity and which of the 
couplings correspond to relevant ones. Most importantly, we cannot draw any strong conclusions from the evidence 
for a fixed point in the case of quantum Einstein gravity |20H28j . As the symmetry changes from full diffeomorphisms 
to transverse diffeomorphisms, and the field content is restricted by the unimodularity requirement, the two theory 
spaces are different. The existence and properties of possible fixed points can therefore be different in the two cases, 
as discussed in ^5] . 
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B. Comments on the equivalence between quantum Einstein gravity and unimodular quantum gravity 

In the following, we will refer to the quantum theory with a full dynamical metric and full diffeomorphism symmetry 
as quantum Einstein gravity, also in cases where the underlying action is not the Einstein-Hilbert action. In the 
unimodular case, = e implies that the volume term is not an operator any more, but simply a number. As such, 
it will be dropped from the action, and no quantum fluctuations contribute to the running of the prefactor, i.e., the 
cosmological constant. In other words, the cosmological constant is removed from the unimodular theory space. In 
principle, one could also attempt to impose the unimodularity condition employing a Lagrange multiplier, in which 
case the cosmological constant would remain a running coupling, presumably resulting in an inequivalent quantum 
model [5]. In quantum Einstein gravity, one can impose = e as a gauge. This does of course not remove the 
cosmological constant from the theory space. This becomes particularly important when the cosmological constant 
corresponds to a relevant coupling at an RG fixed point. Then the predicted values of all other couplings in the infrared 
will depend on the corresponding free parameter, associated to the cosmological constant. On the other hand, in the 
unimodular case, the value of the cosmological constant, which will first enter the theory at the level of the equations 
of motion, is also a free parameter. However none of the predictable values of the irrelevant couplings depend on that 
parameter. Furthermore, as pointed out in |29j . choosing = e as a gauge in quantum Einstein gravity implies 
that a corresponding Faddeev-Popov ghost contributes to the Renormalization Group flow. This additional ghost is 
completely absent in unimodular gravity. Imposing unimodularity directly on the allowed configurations of the metric 
in the path-integral also changes the symmetry from diffeomorphisms to transverse diffeomorphisms; again implying 
differences at the quantum level. Choosing a classically equivalent formulation of unimodularity can allow to keep 
full diffeomorphism invariance [7], however we focus on the other case here. 

Finally, imposing = e in unimodular gravity also implies that the spectrum of quantum fluctuations, i.e., the 
off-shell part of the propagator of metric fluctuations, differs. We will again see this explicitly in Sec. |III C 2[ By 
itself, this already changes the RG flow. We conclude that an equivalence between unimodular quantum gravity and 
quantum Einstein gravity is not be expected. 


C. Deriving the flow equations 

In this section, we will detail the derivation of the RG flow equation for unimodular quantum gravity. We will focus 
on a setting where unimodularity is implemented as a restriction on the allowed configurations in the path integral. 
This has the advantage that it reduces the number of fluctuating gauge degrees of freedom in the path integral, and 
could therefore be expected to yield better results already in simple approximations of the full path-integral. 


1. Wetterich equation 

To examine whether asymptotic safety is realized in unimodular gravity, knowledge of the non- perturbative beta 
functions is required. Here, we will use the framework of the functional Renormalization Group to obtain these. 
In that setting, an infrared cutoff function Rk[—T>^), (with T> denoting a placeholder for the appropriate covariant 
Laplacian) depending on the momentum scale k, is included in the generating functional, which suppresses quantum 
fluctuations with covariant momenta —T)^ < [301 El]. For instance, a theta-cutoff takes the form [33| 

= (/,2 _ (_ 252 )) 0 (fc 2 _ ( 9 ) 

This allows us to obtain the scale dependent effective action T k , which encodes the effect of high-momentum quantum 
fluctuations. Its scale dependence is governed by the Wetterich-equation 

dtVk = isTr (rf + Rk) dtRk, dt = kdk- (10) 

( 2 ) 

Herein T), is the second functional derivative with respect to the quantum fields, which is matrix-valued in field 
space. The supertrace implies a summation in field space, including a negative sign for Grassmann-valued fields. 

/ n'\ 

It also encodes a summation/integration over the discrete/continuous eigenvalues of the kinetic operator T), . For 
reviews, see, e.g., [31]. In the case of gravity, M. Reuter has pioneered the application of the Wetterich equation in 
[25j . Gravity-specific reviews of the asymptotic safety scenario and the application of the functional Renormalization 
Group can be found in, e.g., [35] . 

The application of Renormalization Group methods requires us to set a scale, which seems a challenging task in 
quantum gravity, where no background spacetime is assumed to exist. Here, the background field method |36| can 
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be employed to circumvent this problem: Splitting the full metric into a background and a fluctuation piece provides 
a background covariant derivative. This can be used to define ”high-momentum” and ’’low-momentum” quantum 
fluctuations by decomposing the quantum field into eigenfunctions of the background Laplacian and sorting them 
according to their eigenvalue. At the same time, admitting fluctuations of arbitrarily large amplitude implies that we 
can still perform the functional integral over all metric configurations, as long as the topology is kept fixed. While 
quantum Einstein gravity admits a linear split of the metric, unimodularity requires us to use a non-linear split of 
the form 


9fiu — T “f “t” — 9fii' T -t- ..., (H) 

where the background metric is used to lower and raise indices at each order of the expansion in the fluctuation 
field hfiv We then take the path-integral over the fluctuation field as the definition of the generating functional 
for quantum gravity. (The effect of using such a decomposition in quantum Einstein gravity has been studied in in 
[371135] and [33].) Note that a similar decomposition has to be invoked for calculations involving fluctuations, e.g., in 
a cosmological setting, in the context of (semi-) classical unimodular gravity. 

We then have at our disposal the background covariant Laplacian —D^ which we can use to set up a regulator 

confused with the Riemann tensor) for the fluctuation field 

+ aR)Kx. ( 12 ) 

where a is the prefactor of a possible additional dependence on the background curvature R. While the action is 
invariant under a simultaneous transformation of the background and fluctuation piece {gp.K. —t 5 /iAexp (7 )^, —t 
+ 1/2[/i.,,7..]a.k - l/12[/i.,, + l/6[7.., + ■•■): tlie regulator term is clearly not. The same 

will be true for the gauge-fixing, as we will employ a background gauge fixing here. This implies that background- 
field couplings and fluctuation-field couplings will satisfy different flow equations. Eor instance, the RG flow of the 
background Newton coupling will differ from that of the prefactor of the term quadratic in the fluctuation field and 
in derivatives. In this work, we will not resolve this difference, but instead identify background and fluctuation-field 
couplings, and leave the next step to future work. Thus, while a full calculation would feature beta functions for the 
background and fluctuation couplings where the nontrivial terms can only depend on the fluctuation couplings, our 
approximation will involve a nontrivial dependence on the background couplings. 


2. Second variation of the action 


Our truncation is given by 


Tk = J d'^xefiR). (13) 

Within quantum Einstein gravity, an analogous truncation has been considered in |20H24j . 

The unimodularity condition implies that the number of possible terms in the variation of Eq. ( |13| ) will be reduced, 
since terms which are present in quantum Einstein gravity, such as {6^y/g)f{R) and {5^)f'{R)5R, do not exist here. 
In fact 


6^Tk = j d^xe {f{R)S^R + r{R){SRf) . 


(14) 


Evaluating the variation, starting from the relation Eq. 0, and using a 4-sphere for the background field configu¬ 
ration^, we obtain 


= - / d^xe 


f”{R)h^,D^^D''D-D^Kx + f'iR) ( -^RKuh^'' - ^ 


(15) 


Herein, i.e., we employ a single-metric approximation from now on. 


2 For this configuration, we have f and ^ - S^iA9ck). 
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As a next step, we insert a York-decomposition of the fluctuation field into a transverse traceless tensor, a transverse 
vector, and a scalar (corresponding to the longitudinal vector mode). Note that in contrast to the usual case, there 
is no trace mode, i.e.. 




(16) 


where = 0, = 0 and = 0. 

It turns out that the second variation evaluated on the transverse vector mode vanishes. In other words, the 
dynamics of the vector mode is arising from the gauge-fixing term only, i.e., it is “pure gauge”. This is another major 
difference to the case of f{R) truncations in quantum Einstein gravity, see, e.g., [5DH11]. 

For the transverse traceless tensor mode we obtain 


- / d xehu^r. 


XT 


( 2 ) fiyKX 


h.x = \j d^xef'{R)hl'^ (y-D^ - 


(17) 


Finally, the scalar mode is governed by the following dynamics 


- / d'^xeaT^^la = - d'^xea 


f'{R) ( —RD'^D^ - —D'^D'^D'^ 
^ ^ 16 16 


+f''{R) ( + ^RD^D^D"^ + ^R^D^D^ 


a. 


(18) 


As usual, no mixed contributions To-i, etc. can exist because of the transversality and tracelessness of and Vfj,. 


3. Gauge-fixing 

We choose a gauge-fixing that is related to the harmonic gauge condition, but modified such that it satisfies 

= 0, (19) 

for the spherical background. Accordingly, this choice of gauge fixing only imposes three instead of four gauge- 
conditions, i.e., it only fixes the transversal diffeomorphisms, infinitesimally defined by 

Sodni' = with = 0. (20) 

Note that in models of gravity which are invariant under transverse diffeomorphism, an additional scalar mode 
appears upon linearization. As noted in m, this mode is absent in two cases: If the symmetry is enhanced to full 
diffeomorphism symmetry, yielding standard Einstein gravity, or if the metric determinant remains fixed. Then the 
additional scalar, which plays the role of the determinant, is removed from the model. 

Gauge-fixing only the transverse diffeomorphisms is achieved by using the longitudinal and transversal projectors 
dehned in sn, which read 


Hlm. = -D^{-D^) (21) 

llT/iiy = gij.1^ — (22) 

As they should, these satisfy = IIl^k, 11^^1,11^^ = 0 and TIt k — uk- We now project the 

harmonic gauge on its transversal part |42j and dehne 

(23) 

It is then straightforward to see that g^’'D^Ff^ = 0. Accordingly only three conditions are imposed on the fluctuation 
field, which one can easily see by inserting the York decomposition: It turns out that the gauge fixing does not impose 
a condition on cr, but only on v^, which has only three independent components. These turn out to be gauge modes. 
Indeed 

Fu = ^ [d^ + f) V,. 


(24) 
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Thus the gauge-fixing action reads 


^gfy ~ 


a 






Finally, the Faddeev-Popov ghost action is obtained in the usual way and reads 


Sgh = — I d'^xe i ^ 


j 


(25) 


(26) 


where we have already identified and = 0 = (As we only evaluate the ghost loop contribution 

to the running in the gravitational background couplings, this is already allowed at this stage.) 


4- Jacobian and auxiliary fields 


The York decomposition implies the existence of a nontrivial Jacobian in the generating functional |27j . Here, we 
will deal with this Jacobian by employing the following strategy: From the structure of Sgfy it is obvious that a part 
of the Jacobian can be cancelled by the field redefinition 





(27) 


Employing this field redefinition results in Sgfy = —1/a f d'^xev^ i.e., the vector mode does not 

contribute to the RG flow if we impose Landau gauge. In principle, we could choose to nevertheless impose a 
regulator on that mode with a dependence on the gauge parameter. Here, we take the point of view that a vanishing 
(unregularized) propagator allows us to trivially integrate out the v mode in the path-integral, such that it does not 
affect the effective action. (Alternatively, a gauge-choice of the form = 0 could also be imposed, as in [55], also 
leading to a vanishing contribution of the vector mode.) 

On the other hand, a corresponding redefinition of cr in order to absorb the remaining part of the Jacobian would 
not lead to a simple form of the inverse propagator. Accordingly we introduce auxiliary fields to take into account 
that part of the Jacobian. The corresponding action is given by 


Tk 


aux 





(28) 


where y is a complex Grassmann field and ^ is a real scalar field. These give the same contribution to the flow 
equation, with a relative factor of —2. Accordingly, the total contribution comes with a factor —1/2. 


5. Choice of two regularization schemes and evaluation of traces 


We will study two different regulators in the following. For the first choice, we follow [^ and employ regulators 
which essentially substitute the following covariant Laplace-type operators by k^: 

2R R A A ^ 

A 2 = A + —, Ai = A--, Ao = A--, (29) 

for the transverse traceless tensor, the transverse vector, and the scalar. Herein —D^ = A. 

We choose a Litim-type cutoff [33] 

RkTTi = - ^2)0{k^ - A 2 )), (30) 

for the transverse traceless tensor. Note that the negative sign is exactly as it should be, as in the simplest case 
f{R) = R- For the scalar mode we obtain the slightly lengthier expression 


Rk a 1 


f'iR) 


- Ag) + ^Rik^ - Ag) + _ ^4)^ 


+nR) 


e{e - Ao), 


(31) 








such that for — Aq > 0 


.(2) 


+ Rkal f"{R) 


16 


16 


+ + —fc® + f{R) —i?2 L' 


48 


16 


(32) 


Since the function f{R) appears explicitly, its scale-derivatives will feature on the right-hand side of the Wetterich 

equation. They will lead to a schematic structure of the form dtg = cig^ + C 2 dtg -I- ... of the flow equation for the 

2 

couplings. This results in nonperturbative resummation structures, i.e., dtg ~ i^c 2 ' prefactor C 2 can contain 

further couplings, this choice of ’’spectrally adjusted” [3311111 regulator roughly corresponds to a resummation of an 
infinite series of polynomial terms in the couplings. 

Explicitly, we will use the following derivatives 

dtf{R) = ( 2 / + dtf - 2i?/") , 

dtf'iR) = dtf -2Rf''. 

in terms of the dimensionless function /()?) = k~‘^f{R), where R= 

For the auxiliary fields and ghost we take 

i?fcghi = -y2(/c"-Ai)0(fc2-Ai), (35) 

i?fcauxi = ^ - Ag + I - Ao)) e{e - Ao). (36) 


(33) 

(34) 


In order to test the reliability of our results, we will actually perform a fixed-point search with two different regu¬ 
larization schemes. As our second choice we employ a regulator, which essentially substitutes covariant Laplacians 
A = —D^ with k^ in the regularized propagator. The additional curvature dependence which is introduced in the 
regulator when using the operators Eq. (29) is removed in this choice, resulting in a shift in possible poles of the flow 
equation. With A = —D'^, this choice corresponds to 


RkTT2 = -lfiR)ik^ - A)0(fc2 - A)), 


(37) 


for the transverse traceless tensor. For the scalar mode, we choose a regulator of the form 


Rk a, 2 


'f{R) - A^) + ^(fc® - A®)^ 

+nR) - A") - - A®) -f ^(fc4 


6l(fc2 


A). 


For the ghost, we choose 


(38) 


Rkgh,2 = V2i-k‘^ + A)e{k‘^-A). (39) 

The auxiliary fields come with a regulator of the form 

i?fcaux,2=-^ (^-fc4 + A2 + |(fc2-A)^0(fc"-A). (40) 

In both regularization schemes, we sum over the eigenvalues of the corresponding Laplacians, which can be obtained 
for both A and A^ from the following table, where the multiplicities are not affected by the curvature-dependent shift 
between A and A,,. 

To convert the sum over eigenvalues into an integral, we employ the Euler-MacLaurin formula. In this step, 
additional terms which depend on the derivatives of the integrand at the lower boundary, arise. No contributions 
from the upper boundary exist, as 9{k'^ — x) — 0 for x —>■ 00 . From the lower boundary, only the first few terms 
contribute: As dtRkff + Rk) is a polynomial of finite order in the eigenvalues, only the lowest few orders in 
derivatives, when evaluated at the lower boundary, can contribute. 
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eigenvalue 

multiplicity 


n(n + 3)-4 p „ , 

(n + 2)!(2n+3) 


12 ••• 

6n! 

A, 

n(n + 3)-4 p . ^ 

(n+l)!n(n+3)(2n+3) 


12 -^5 -^-5 ••• 

2(n + l)! 

Ac. 

n{n+3) p 0 Q 

5(n+l)!(n+4)(n—1) (2n + 3) 



6(n + l)! 


TABLE I: Eigenvalues and multiplicities of Laplacians acting on transverse traceless tensors, transverse vectors, and scalars, 

m- 


IV. RESULTS: ASYMPTOTIC SAFETY IN f{R) 
A. Flow equations 


The flow equation for f{R) reads 

dtf{R) = -4/(f?) + 2i?/'(f?) + ^^(J-TT+J'.+J'gh + .Faux), (41) 

where the contributions from transverse traceless tensors RtTi scalars auxiliary helds J^aux and Faddeev-Popov 
ghosts Jgh depend on the choice of regulator, as discussed above. 

For the first choice of cutoff, which is a function of A^, we obtain 


^tt — 77 : 


89 60 40 


1 / 89 20 20 127i? 


Ra = 


18 i ?2 Jij 2 /' y 18 i ?2 Ji 567 

1 


( 2 /' + dtf - 2 /"i?) , 




+ ( 2 /' + dtf - 2 /"i?) 

+ (aj'" - 2/"'7?) 
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(42) 


(43) 

(44) 

(45) 


As a main structural difference to quantum Einstein gravity with a full dynamical metric, including a conformal 
factor [20H22] . one should note that the function / does not appear on the right-hand side of the flow equation, i.e., 
the cosmological constant is not part of the theory space. 
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For the second regulator, that imposes cutoffs on A, we obtain 
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Ftt = 
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90i?2(-3 + i?) 
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(47) 

(48) 

(49) 


Comparing the two equations for dtf{R), we note that the different choice of regularization scheme mainly serves 
the purpose of changing the singularity structure of the equation, which can play a role in the search for global 
solutions |22j . Apart from that, the structure is similar in both cases, with differences only in the numerical prefactors 
of most terms. 


B. Fixed points 


To search for fixed points, we expand f{R) in a polynomial around vanishing curvature, where the cosmological 
constant is again absent. 


10 


/(/) = ^ a„/" 


(50) 


Note that the Newton coupling is given by G = Thus a negative fixed-point value for oi translates into a 

positive microscopic Newton coupling. Using the first regularization choice Eq. (421 - Eq. (451, we find a number of 


fixed points, where we list only the most stable one in Tab. [n| Beyond iR it becomes rather cumbersome to evaluate 
all solutions of the fixed-point equations. We thus pick the most stable fixed point that exists in the four truncations 
up to 04 . At R? and beyond we perform a numerical search for the solution of the fixed-point equation in the vicinity 
of the fixed-point coordinates at lower truncation order, and no longer investigate all solutions. 

The fixed point at n = 1 is related to that found in [^, which neglected the additional terms arising from the 
Euler MacLaurin formula. We clearly observe that the fixed point at truncation order n + 1 is a continuation of the 
fixed point at order n, as the coordinates of the fixed points lie reasonably close to each other. This stability under 
extensions of the truncation indicates that this is not an auxiliary fixed point. Interestingly, the fixed-point values 
of the couplings a„, n > 5 are four orders of magnitude smaller than the highest-order coupling. This suggests that 
an approximation of the fixed-point action employing only the first few terms could be reasonable, for instance if 
cosmological consequences of unimodular asymptotic safety are deduced from ”RG improved” calculations. 

The corresponding critical exponents are given in Tab. m and again seem reasonably stable under extensions of 
the truncation, cf. Eig. Similarly to the case of quantum Einstein gravity, two exponents with positive real part 
are observed, i.e., R and R^ form relevant directions. In this approximation, the two critical exponents with positive 
real part form a complex conjugate pair, which could hint at a necessary extension of the truncation to obtain better 
numerical precision. Higher-order operators become increasingly irrelevant. The anomalous dimension is a positive 
contribution to the critical exponents for all operators, i.e., quantum fluctuations shift all operators towards relevance. 
Comparing the value of the critical exponents 9i beyond i = 3 to the canonical dimension di = — (2i — 4), we observe 
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TABLE II: Fixed-point values for the dimensionless couplings a„ in a polynomial expansion around vanishing curvature, at 
increasing truncation order. 
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2.122 ± i 1.232 


2.778 ± i 1.232 

-1.233 


2.832 ± i 0.781 

-1.113 

-3.111 


2.912 ± i 0.687 

-1.172 

-3.415 

-5.235 


2.863 ± i 0.654 

-1.155 

-3.328 

-5.447 

-6.997 


2.862 ± i 0.671 

-1.197 

-3.382 

-5.380 

-7.464 

-8.863 


2.847 ± i 0.673 

-1.204 

-3.418 

-5.483 

-7.325 

-9.398 

-10.708 


2.841± i 0.678 

-1.221 

-3.440 

-5.530 

-7.534 

-9.236 

-11.324 

-12.698 


2.834 ± i 0.681 

-1.228 

-3.462 

-5.569 

-7.585 

-9.545 

-11.093 

-13.224 

-14.514 


TABLE III: Critical exponents at the fixed point shown in Tab. [11] The first line is of course understood to feature only one 
real critical exponent. 


a decreasing distance up to 6q, cf. Fig. Beyond, our truncation is not large enough to produce convergent results, 
as can be inferred from the comparison between the largest and second-largest truncation. Tentatively extrapolating 
this trend would indicate that canonical scaling could be recovered at large powers of the curvature. 

We observe an interesting difference to quantum Einstein gravity, where the gap between the smallest relevant and 
largest irrelevant eigenvalue is Aj/v = 5.58 [24]. In our case, we find Af/yunimod ~ 4.06, i.e., a significantly reduced 
gap. This is related to the fact that canonically irrelevant operators are shifted further into irrelevance in quantum 
Einstein gravity, whereas they are shifted into relevance in unimodular asymptotic safety. 

At this stage, we should comment on the approximations we have used. There are several sources of uncertainty in 
our results: Firstly, using a spherical background to project on the function /(i?) entails ambiguities. Operators with 
the same number of derivatives are projected onto simultaneously: On a sphere, derivatives of the curvature tensor 
and its contractions vanish. Further, the Weyl tensor vanishes, allowing to re-express the Riemann tensor in terms of 
the Ricci tensor and the Ricci scalar. Thus, the remaining ambiguity lies in the fact that on a 4-sphere 
Accordingly we derive the spectrum of fluctuations, i.e., F^^^ from an action of the form f{R), without contributions 
from other tensor structures. Using a sphere to evaluate the traces and project onto powers of R implies that we 
project on R^ + ciR^i^R^'^ at fourth order in derivatives, and so on, with unknown coefficients cy This ambiguity 
is inherent in the setup of our calculation and is owing to a compromise between uniqueness of the projection and 
computability. Here, our calculation is on a par with similar challenges encountered by f{R) approximations in 
quantum Einstein gravity. 

To address this ambiguity, we make a useful observation: Had we considered beta-functions for the disentangled 
system which distinguished between R^ and R^^R^'^^ etc., a fixed point in that system necessitates the existence of a 
fixed point in our simplified system. Furthermore, a relevant direction in the simplified system can only occur, if (at 
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FIG. 1: Left panel: Real part of the five largest negative criitical exponents (the first two are degenerate) and canonical 
dimensions of the couplings (dashed lines). Right panel: Difference between the critical exponents 6 i (starting with the third- 
largest) to the canonical dimension di = —{2i — 4) for the truncation including np to Iv (dark red dots) and up to (blue 
dots). Beyond the dashed purple line our results have not converged. 
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TABLE IV: Fixed-point values using the second regularization scheme. 


least) one relevant direction exists in the disentangled system. Thus, our estimate of the number of relevant directions 
provides a lower bound on the number of relevant directions in the disentangled system, and a discovery of a fixed 
point in our approximation is a necessary condition for the existence of a fixed point in the disentangled system. 

Secondly, we have truncated the theory space very severly. This is important in two places: Firstly, our truncation 
is not closed in the sense that operators beyond those present in the truncation are generated on the right-hand-side of 
the Wetterich equation. Secondly, further operators also change the spectrum of fluctuations, i.e., T^^^ which enters 
the evaluation of the beta-functions for the coefficients at. 

Thirdly, we have employed a so-called single-metric approximation. This consists in neglecting the fact that no 
background coupling should appear on the right-hand side of the Wetterich equation, as the functional derivatives 
necessary to obtain T^^^ are taken with respect to the fluctuation field. Thus, fluctuation field couplings should 
appear on the right-hand side of the Wetterich equation. Performing the necessary evaluation of fluctuation-field and 
background-field flows at the level of an f{R) truncation is beyond the scope of this paper. 

All these approximations can affect the stability of our results. One could expect that if the errors introduced into 
the results by the approximations are small, then the regulator-dependence should not be too large. In this spirit we 
study whether a fixed point with similar properties can also be found using the second regularization scheme, Eq. (46 )- 
Eq. (49). We find a fixed point with similar properties (see Tab. IV and [v|), thus strengthening the conclusion that 
the fixed point is not a truncation artifact. 

The fixed point found at positive Newton coupling with the second regularization scheme shares several charac¬ 
teristics with the fixed point found in the other scheme. Firstly, there are two positive critical exponents which are 
comparable to the real part of the two relevant critical exponents in the other scheme. We observe that fixed-point 
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2.441 

11.349 

-1.740 


4.683 

1.498 

-1.952 

-3.444 


4.542 

1.353 

-2.099 

-3.611 

-5.329 


3.937 

1.006 

-2.460 

-3.708 

-5.076 

-7.156 


3.899 

1.204 

-2.646 

-3.656 

-5.209 

-6.923 

-9.058 


3.969 

1.296 

-2.653 

-3.642 

-5.414 

-7.051 

-8.854 

-10.999 


3.953 

1.293 

-2.711 

-3.647 

-5.388 

-7.224 

-8.936 

-10.791 

-12.947 


3.933 

1.289 

-2.762 

-3.653 

-5.385 

-7.192 

-9.089 

-10.851 

-12.746 

-15.907 


TABLE V: Critical exponents at the fixed point shown in Tab. |IV[ 


values and critical exponents converge slower than using the first regularization scheme, which could suggest that 
the first scheme is actually better adapted to the properties of unimodular f{R) gravity. Secondly, the fixed-point 
values for the couplings a„, n > 4, are at least four orders of magnitude smaller than the leading fixed-point values. 
Comparing the fixed-point coordinates in the two regularization schemes, one should keep in mind that all but the 02 
coupling have nonvanishing canonical mass-dimensionalities, i.e., the fixed-point values of all couplings should not be 
expected to be universal. 

With only two relevant directions, the number of free parameters is significantly smaller than that implied by 
the analysis in |46j . with a considerable smaller truncation, where unimodularity was implemented together will full 
diffeomorphism symmetry, necessitating the introduction of Stueckelberg fields. This implies a significantly higher 
predictive power of our implementation, where = const is implemented as a restriction on the path-integral 
measure. It also clearly shows that while the two implementations are equivalent classically, one should not expect 
them to be so on the quantum level, see also [5]. 

It is now interesting to examine the system while dropping the ” RG-improvement” terms dtf on the right-hand 
side, which arise as we chose a regulator that depends on /. Neglecting these terms can have no effect on the positions 
of the fixed points, but can change the value of the critical exponents. We observe that the signs of the critical 
exponents are stable - another indication for the robustness of our result - but their numerical values can change 
quite a bit. This is actually similar to the results in quantum Einstein gravity: Using the flow equations derived in 
[101122], but dropping the RG improvement terms can produce different values for the critical exponents. 

We conclude that the anomalous dimensions, related to dtf, are rather large. This already suggests that it will be 
important to examine the system without employing a single-metric approximation, as this will allow to disentangle the 
background-beta functions from the anomalous dimensions, which are related to the fluctuation field^. As discussed 
in [47], a resolution of the difference between background field and fluctuation field could also play a major role in 
obtaining viable global solutions. 

Our work also has implications for the asymptotic safety scenario in full gravity: It clarifies that the existence of 
a fixed point with three relevant directions in f{R) truncations is not due to the dynamics of the conformal mode. 
While we in fact find one relevant direction less in unimodular gravity, this is due to the fact that we work in a 
different theory space, arising from a different definition of the degrees of freedom and the symmetries. In fact, 
simply dropping the contribution from the conformal mode in full gravity will not result in the removal of the relevant 
direction resulting from the presence of the cosmological constant. It is the definition of a different theory space which 
does not contain the cosmological constant, which lowers the number of relevant directions. Further, our result can be 
read as a further confirmation of asymptotic safety in quantum Einstein gravity: If the removal of only the conformal 
mode could destabilize the existence of a fixed point, one might conclude that the fixed point in quantum Einstein 
gravity was just a truncation artifact. 


^ For a bimetric study of the gravitational Renormalization Group flow that uses the decomposition Eq. that is also required for the 
unimodular case, see m- 
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C. Global solutions 


To analyze whether the fixed-point equation for f{R) might admit global solutions, we bring it into normal form, 
i.e., we rewrite it such that the highest derivative, f'", occurs with a prefactor normalized to one. It is a third-order 
differential equation, and thus a solution comes with three free parameters. Singularities in the equation decrease the 
number of parameters, as they impose constraints on the solution, if it is to be continued through the singularity. 
Let us first consider the second choice of regularization scheme, imposing a cutoff on A. In that case, inspection of 
Eq. (46)- Eq. (49) reveals singularities at i? = 0, i? = 4 and R = — 6 . Additional singularities arise from the prefactor 
of the term ~ /"', which has real zeros at i? « —0.916 and R « 1.624. A seeming singularity at 7? = 3 is actually not 
there, as it also occurs in the denominator of the prefactor of the term ~ f'". So it does not in fact show up when we 
bring the equation into normal form. Note that while the equation for /(i?) was derived using a spherical background, 
we postulate its validity for R < 0 here, and include singularities at negative R in the parameter counting. Taken 
together, these singularities already overconstrain the system, such that a global solution is unlikely to exist. This 
situation is very similar to that found in |23) analyzing the equation derived in | 2 ()j . As we have essentially followed 
a very similar procedure here, and in particular imposed the cutoff in a similar way, it is not unexpected that we find 
a differential equation with similar behavior. Following the insight in |22) . we can derive an alternative equation by 
adapting the choice of regulator, and imposing a cutoff on combinations such as A -|- ^, instead of on A, which in 


fact corresponds to our first choice of regularization scheme. In that case, we should consider Eq. (42|-Eq. (45) to 


determine the existence of fixed singularities. We again observe a fixed singularity at i? = 0. Additional singularities 
again arise from the prefactor of the /"' term, and lie at i? « —0.56 and R « 1.55. At a first glance, another 
singularity seems to lie at i? = —3, but again does not exist in the normal form of the equation. We conclude that this 
third-order equation exhibits precisely three singularities, which suggests that a global solution could exist. Whether 
the asymptotic structure at large R imposes further constraints, requires a more detailed analysis which we leave to 
the future. Here, we only observe that an asymptotic solution is given by f{R) ~ AR^ at leading order, as expected 
from canonical dimensionality. 


V. TOWARD THE REAL WORLD: ADDING MATTER 

If a theory of quantum gravity is to be applied to our universe, dynamical matter degrees of freedom must be 
accounted for. Here we take the point of view that matter fields must be included at the microscopic level. Thus, 
our truncation should not only include dynamics for gravity, but also fermions, vector bosons, and scalars. We will 
consider only minimally coupled matter, and disregard further interactions for the moment. This is an approximation, 
as metric fluctuations induce matter self interactions in the UV [48l 09]. Interestingly we can take advantage of 
structural similarities between quantum Einstein gravity and unimodular gravity: The matter loop contribution 
to the background Newton coupling beta function is the same in the two settings. This will change within more 
sophisticated truncations. For instance, it is no longer true when wave-function renormalizations for the matter fields 
are included, as these are derived from different graviton-matter-vertices. We can thus take the matter contributions 
from [SOI mi and add these to the beta function for the background Newton coupling, where we use the n = 1 
truncation and identify ai = —l/( 167 rG) and subsequently expand to second order in G. We then have that 

Pg = ‘2G-^ (20 - 2Nd -Ns+ mv + Mrs) , (51) 

bTT 

where Ny is the number of abelian vector fields. Nr the number of Dirac fields. Ns the number of real scalars and 
Nrs the number of Rarita-Schwinger fields with spin 3/2. Herein a type-H cutoff (nomenclature as in [3D]) has been 
used for the matter fields, as it is required for the proper treatment of fermions [52] . For the gravitational contribution, 
we have used the cutoff imposed on A 5 . In analogy to the Einstein gravity case, there is a bound on the number of 
Dirac fermions and scalars, if the number of vectors and spin 3/2 fields is fixed. Disregarding Rarita-Schwinger fields, 

2Nd+Ns <20 + mv, (52) 

is required for a positive Newton coupling at the fixed point. While the microscopic value of this coupling is not 
(yet) restricted observationally, its infrared value must be positive. In all known truncations, the sign of the Newton 
coupling is preserved under the RG flow. We therefore exclude negative fixed-point values based on the positivity of 
the observed Newton coupling in the infrared. 

Considering supersymmetric models, pure supergravity with a gravitino admits a viable fixed point. On the other 
hand, the matter content of the minimal supersymmetric standard model {Ns = 49, Nr = 61/2 and Ny = 12) seems 
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FIG. 2: We show the bound on the allowed number of Dirac fermions No and scalars Ns with Nv = 12 and Nas = 0. The 
blue thick line is obtained with the first regulator, and corresponds to Eq. (511, while the thin dashed blue line corresponds to 
the second regulator. For comparison, the purple dashed line is the bound in quantum Einstein gravity from m- The green 
dot corresponds to the Standard Model degrees of freedom, and the turqoise asterisk additionally contains three right-handed 
fermions (neutrinos) and a dark matter scalar. The red dot, corresponding to the degrees of freedom of the MSSM lies outside 
the allowed region. 


excluded. Whether this result will be affected when the truncation includes supersymmetric interactions, and also an 
appropriate regularization scheme [53] is used, is presently unclear. 

Most importantly, the Standard Model matter content {Ny = 12, Ns = 4, Njj = 45/2) admits a gravitational 
fixed point at G* > 0. We can also add several additional scalars or fermions, that could constitute dark matter 
and still find a viable fixed point. This result is confirmed if we use a regulator imposed on A = —D^ for the 
gravitational contribution. This changes the factor 20 in Eq. (51) to 14 and yields a tighter bound on fermions and 
scalars. Crucially, the bound still includes the Standard Model, cf. Fig. This indicates that unimodular quantum 
gravity can be compatible with the matter content of our universe. 

While this result is obtained within the simplest truncation, and requires much further investigation, it is encour¬ 
aging that we cannot find a way to exclude our model at this stage from phenomenological considerations of matter. 
On the other hand, we actually observe a tighter bound on matter than within quantum Einstein gravity. This arises, 
as the bound comes from a balance of gravitational and vector with fermion and scalar fluctuations. As the gravita¬ 
tional contribution is smaller in the unimodular case, the balance is reached for lower numbers of fermions/scalars. If 
this result persists beyond our truncation, future discoveries of additional (dark matter) fermions and scalars might 
rule out unimodular asymptotic safety (or necessitate the introduction of additional vectors with the corresponding 
scalar modes to make them massive). Of course, much more detailed studies are necessary in order to investigate 
the reliability of these bounds, and also understand inhowfar these are universal. Crucially, the form of the bound 
presented in this paper should be understood as coming with considerable systematic errors. Beyond its implications 
for unimodular quantum gravity, the bound Eq. (52) nevertheless points towards a potential option how experimental 
constraints might be put on quantum gravity models from low energy observations. 


VI. CONCLUSIONS 

We have studied the Renormalization Group flow of unimodular quantum gravity and found an interacting fixed 
point, further strengthening the evidence for unimodular asymptotic safety. 

In particular, we have used a truncation of the effective action to a function of the curvature scalar, f{R), and 
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derived the flow equation for this function. We have subsequently expanded in powers of the curvature up to and 
found a fixed point with two relevant directions, i.e., two out of ten couplings correspond to free parameters of the 
model. 

We have studied two different regularization schemes, and found a fixed point with similar properties and two 
relevant directions in both schemes. While this provides evidence that the fixed point is not an artifact of the 
approximation scheme, the remaining scheme dependence indicates that further studies are necessary to determine 
the fixed-point properties with quantitative precision. 

As a main difference to quantum Einstein gravity, the cosmological constant is not part of the theory space in the 
unimodular setting. The corresponding fine-tuning problem is therefore avoided. This makes unimodular asymptotic 
safety attractive from a phenomenological point of view. The cosmological constant will only enter, once the RG flow 
has been integrated into the infrared, where the effective equations of motion can be derived. As we have shown, 
these can then be reformulated with the help of the Bianchi-identities, and the cosmological constant appears as a 
constant of integration. 

Furthermore the unimodularity condition = const reduces the number of propagating components of the metric 
in the quantum theory and most importantly removes the conformal mode from the path integral. This avoids the 
corresponding instability of the Euclidean path integral in an Einstein-Hilbert truncation. Besides, the reduction in 
degrees of freedom (by which we here mean components of the quantum field and not necessarily physical degrees of 
freedom) results in computational simplifications. Most importantly, we find that an appropriate gauge fixing removes 
all but the transverse traceless and one scalar mode of the metric. 

Performing a first step towards phenomenology, we have also studied the compatibility of unimodular asymptotic 
safety with minimally coupled matter degrees of freedom. Here we make use of structural similarities between quantum 
Einstein gravity and unimodular gravity at the level of the background couplings, which imply that the matter 
contributions to the running of the background Newton coupling in unimodular gravity agree with [50115T] . We find 
bounds on the number of allowed scalars and fermions, at a fixed number of Abelian vector bosons. At this order 
of the approximation it is only the gravitational contribution to the Renormalization Group flow that differs from 
the result in quantum Einstein gravity, and which yields a slightly tighter bound on the allowed number of matter 
fields. Most importantly, we find that within this approximation, the Standard Model, as well as small extensions by, 
e.g., a dark matter scalar, can be accommodated in unimodular asymptotic safety. This suggests that unimodular 
asymptotic safety could pass an important observational test. Moreover, this provides a first example, where different 
models of quantum gravity might be distinguished observationally: As - within simple truncations - unimodular 
asymptotic safety and quantum Einstein gravity impose different bounds on the number of matter fields, detection 
of certain BSM-matter models, e.g., at the LHG, could potentially rule out one of these models while admitting the 
other. To further quantify this exciting possibility, nonminimal matter-gravity couplings need to be included in the 
truncations. 

From here, several interesting routes are open toward future work: As has recently been emphasized in I371I3H], 
the exponential parameterization of the metric in terms of the fluctuation field, that we employed here, can also 
be advantageous in quantum Einstein gravity. From our result, it is only a small step towards an analysis of the 
f{R) truncation in that parameterization, as only the conformal mode needs to be added and the contribution from 
fluctuations of y/g needs to be added to the other propagators of the other modes. 

Further, an extension of the truncation in a bimetric direction is possible as in quantum Einstein gravity |28j and 
seems indicated, see, e.g., m- This is possible both within pure gravity as well as including matter. In that case, 
it is particularly interesting to investigate whether the unimodularity condition on the matter-graviton vertices can 
have observable consequences. For instance, one could imagine that differences in the vertices result in a change of 
the number of relevant directions in the gravity-matter sector. This could even lead to differences in predictions for 
the low-energy theory, e.g., along the lines of [51] . 
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